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N*\^)QjC*  ABSTRACT 

The  existing  analytical  solutions  for  gas  lubricated  stepped 
sector  thrust  bearings  are  reviewed.   An  experiment  performed  at  high 
bearing  numbers  (123  to  165)  using  air  at  atmospheric  pressure  is 
described.   Results  are  presented  of  tests  conducted  using  a  bearing 
optimized  by  the  methods  of  Ausman  and  Wildmann  for  high  speed  operation, 
Most  of  the  design  optimizations  are  confirmed  for  high  bearing  numbers; 
the  exception  is  load  carrying  capacity.   The  load  capacity  curve  ob- 
tained takes  the  same  form  as  that  predicted  by  theory  but  is  lower  for 
any  given  bearing  number  by  a  constant  percentage.  This  points  to  a 
need  for  re-evaluation  of  determination  of  load  carrying  capacity. 
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NOMENCLATURE 

2 

A  Bearing  area  (in.  ) 

b  Radius  ratio  ■  r./R 

h  Lubricating  film  thickness  (in.) 

h.  Film  thickness  in  leading  edge  region  (in.) 

h_  Film  thickness  in  trailing  edge  region  (in.) 

L  Load  parameter  ■  W/pa  A 

n  Number  of  sectors 

N  Rotation  rate  (rpm) 

p  Pressure  (psi) 

q  Compression  ratio  -  ni/h9  =  (  A  +  h  )/h-  =  1  +  l/°< 

Q  Mass  flow  rate  (lbm./sec. -in.) 

r  Radial  coordinate  (in.) 

rt  Bearing  inner  radius  (in.) 

R  Bearing  outer  radius  (in.) 

t  Time  coordinate  (sec.) 

V  Velocity  of  fluid  (in. /sec.) 

W  Bearing  load  (lbf.) 

z  Axial  coordinate  (in.) 

°<  Minimum  clearance  ratio  ■  h_/& 

A  step  depth  (in.) 

S  Angular  width  of  groove  between  thrust  pads 
(degrees  or  rad.) 

A  Bearing  number  ■  —  ,?- — 

/\_  Bearing  number  based  on  step  depth  ■    p  a* — ' 

/{  Modified  bearing  number  ■  .A.(|-b*) 

/\_  Modified  bearing  number  based  on  step  depth  =  A-A0""b  ) 


•  Sector  arc  length  (degrees  or  radians) 

If  Step  ratio  -  r/y 

P  Lubricant  density  (lbm./in.  ) 

0  Coordinate  angle  (radians) 

2 

q  Viscosity  (reyns  ■  lbf . -sec. /in.  ) 

f  Arc  length  of  leading  edge  region  (degrees  or  radians) 

^  Angular  velocity  (rad./sec.) 

Subscripts 

1  Refers  to  leading  edge  region 

2  Refers  to  trailing  edge  region 
a  Refers  to  ambient  condition 

r  Radial 

©  Circumferential 

A  Refers  to  the  step 

opt  Optimized 
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1.    Introduction. 

In  a  rapidly  expanding  technology,  which  creates  demands  for  machin- 
ery components  that  will  operate  in  progressively  more  severe  or  unusual 
environments,  it  is  soon  found  that  "conventional"  lubricants  for  moving 
or  rotating  parts  will  no  longer  serve.   Because  they  possess  several 
unique  advantages,  gases  are  a  logical  consideration  as  a  lubricant  when 
seeking  solutions  to  many  design  problems. 

Gas  lubricants  are  particularly  attractive  at  extremely  high  or  at 
very  low  temperatures.   Since,  in  contrast  to  liquids,  the  viscosity  of 
gases  is  relatively  constant  and  even  increases  slightly  with  temperature, 
gas  lubricants  have  an  obvious  application  in  very  high  temperature  rotat- 
ing machinery.   At  the  other  end  of  the  scale,  elements  such  as  nitrogen 
and  helium  retain  their  desirable  lubricating  qualities  at  cryogenic 
temperatures.   Indeed,  the  "boiloff"  from  liquified  gases  is  often  used 
as  the  lubricant  for  expansion  turbine  gas  liquifiers. 

Another  important  advantage  of  gas  lubricated  bearings  is  the  extreme- 
ly low  power  loss  from  friction  at  very  high  speeds.   This,  plus  the  fact 
that  this  power  loss  is  relatively  unvarying,  dictates  the  use  of  gas 
bearings  in  precision  applications  such  as  high  speed  gyroscopes.  The  low 
friction  loss  at  high  speeds  plus  high  temperature  operation  makes  gas 
lubricated  bearings  the  most  desirable  type  for  gas  turbines  and  turbo- 
jet engines. 

In  many  areas  gas  lubricants  can  solve  the  problem  of  contamination. 
When  pumping,  compressing,  or  expanding  high  purity  gases  or  refrigerants, 
the  use  of  the  process  fluid  as  a  lubricant  will  prevent  mixing  of  the 
fluid  with  oils  or  greases.   In  food  processing  machinery,  elimination  of 
petroleum  lubricants  prevents  contamination  of  foods.   In  nuclear  reactors, 
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it  is  desirable  to  eliminate  oils,  greases,  and  solid  lubricants  which 
break  down  or  form  radioactive  isotopes  when  subjected  to  radiation. 

As  with  conventional  bearings,  gas  lubricated  bearings  fall  into  two 
major  classes,  self-acting,  or  hydrodynamic  bearings,  and  pressurized, 
or  hydrostatic  bearings.   There  exists  a  wide  variety  of  both  journal  and 
thrust  type  bearings  in  each  class.   Self-acting  bearings  are,  in  general, 
the  least  complex,  requiring  neither  separate  gas  supplies  and  piping  nor 
fluid  passages  and  orifices  to  complicate  the  design  of  a  machine. 

Unfortunately,  gas  lubricated  bearings  have  only  recently  emerged 
from  the  realm  of  laboratory  curiosities.   Currently  most  types  are  sub- 
ject to  serious  instabilities  and  require  extremely  close  machining 
tolerances  with  very  smooth  finishes.   Air  film  thicknesses  are  so  small 
that  thermal  distortions  are  a  major  cause  of  failures.   To  gain  a  better 
understanding  of  the  capabilities  and  problems  of  gas  lubricated  bearings, 
this  project  was  undertaken. 

Because  of  an  interest  developed  while  the  author  was  working  at 
AiResearch  Manufacturing  Company  on  a  summer  industrial  tour,  this  study 
will  be  confined  to  self-acting  thrust  bearings  of  the  stepped  sector 
type.   There  are  several  arguments  which  support  this  choice. 

Because  the  stepped  sector  thrust  bearing  (See  Figure  1  for  typical 
configuration)  is  a  common,  popular  type  noted  for  its  simplicity  and 
relative  ease  of  manufacture,  a  large  volume  of  definitive  theory  and 
design  data  has  been  developed.   The  most  recent  and  acceptable  theory  as 
yet  lacks  complete  experimental  verification.   The  recent  development, 
at  AiResearch,  of  a  thrust  bearing  test  rig  capable  of  very  high  speeds 
affords  an  opportunity  to  obtain  data  in  regions  of  operation  previously 
unattainable.   Arrangements  were  made  for  the  construction  of  a  second 
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Thrust  Bearing  Test  Rig  ML1606  at  the  Naval  Postgraduate  School  in 
order  to  conduct  this  research. 
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2.   Theoretical  Considerations. 

In  1886  the  English  scientist,  Osborne  Reynolds,  published  the  de- 
fining differential  equation  for  thick  film  lubrication.   Several  years 
later,  Kingsbury  in  the  United  States  and  Michell  in  Australia  demonstrat- 
ed the  applicability  of  the  Reynolds  Equation  to  thrust  bearings  by 
development  of  the  tilting  pad  bearing.   In  1918  Lord  Rayleigh  [1]  found 
that  for  an  infinitely  wide,  incompressible  film,  the  optimum  shape  for 
a  thrust  bearing  is  the  stepped  configuration.   However  the  step  bearing 
can  be  optimized  for  only  one  expected  operating  speed  or  load  but  the 
tilting  pad  bearing  automatically  adjusts  its  pad  angle  to  a  new  optimum 
as  speed  or  load  changes.  Additionally,  when  side  leakage  is  taken  into 
account,  the  stepped  bearing's  load  capacity  markedly  diminishes  and  its 
advantage  over  the  tilting  pad  bearing  disappears.   For  these  reasons 
this  bearing  type  was  little  used  for  liquid  film  lubrication.   It  was 
not  until  1950  that  Archibald  [2]  published  the  first  analysis  of  a 
finite  stepped  thrust  bearing  using  an  incompressible  film. 

When  gas  is  used  as  the  lubricant,  the  step  bearing  becomes  desir- 
able for  two  reasons.   First,  for  a  tilting  pad  bearing  to  retain  its 
continuous  optimum  operation,  the  position  of  the  pad  pivot  point  must  be 
continually  changed.   This  is  the  effect  of  compressibility  on  the  bear- 
ing design.  A  pivot  is  required  that  consists  of  a  cam  surface  that  is 
very  difficult  to  machine.   Secondly  very  close  machining  tolerances  and 
finishes  are  required  on  gas  bearings  and  since  the  tilting  pad  bearing 
is  composed  of  a  number  of  close-fitting  parts,  it  is  extremely  costly 
to  produce.   In  comparison  the  step  bearing  is  relatively  cheap. 

When  the  analysis  takes  into  account  compressible  lubricants,  the 
differential  equation  becomes  nonlinear  and  several  assumptions  become 
necessary  in  order  to  obtain  an  analytic  solution.   Even  though  the 
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resulting  solution  is  approximate  it  does  establish  a  base  for  optimiza- 
tion of  bearing  geometry  and  permits  correlation  of  theory  to  specific 
test  bearings. 

Figure  1  depicts  the  geometry  of  a  typical  stepped  sector  thrust 
pad.   The  upper  thrust  plate  (plane  surface)  rotates  with  constant  angu- 
lar velocity  <^>  while  the  lower  (stepped)  surface  is  stationary.  Radial 
grooves  of  angular  width   8  admit  gas  at  ambient  pressure  to  the  lead- 
ing and  trailing  edges  of  each  pad.   The  pad  itself,  of  sector  arc 
length   V  ,  is  divided  into  two  sectors,  one  sector  of  length  ^  being 
recessed  by  an  amount  a  relative  to  the  other.   The  larger  clearance 
in  the  recessed  region  is  designated  h,   ,  while  the  smaller,  or  mini- 
mum, clearance  is  designated  h2 

In  cylindrical  coordinates,  the  compressible  lubrication  equations 
for  slider  bearings  take  the  form  [3]: 


t-    ~      M   -S-^f-  (1) 


or        •  o  *" 


1  ^p    u    ^-V~  (2) 


^  *  ^-LlS^O  (3a) 

r   or        r    dQ 


(4) 


where  V  and  VQ    are  the  radial  and  tangential  components  of  lubri- 
cant velocity,   Qr  is  the  radial  mass  flow  rate  per  unit  circumferential 
width  and  Q    is  the  circumferential  mass  flow  rate  per  unit  radial  width. 
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To  obtain  these  equations,  the  following  assumptions  were  made: 
(a)  the  flow  is  steady  (does  not  vary  with  time);  (b)  the  film  thickness 
is  at  least  two  orders  of  magnitude  smaller  than  r  or  r©   ;  (c)  the 
axial,  z,  component  of  velocity  and  its  derivatives  are  negligibly  small; 
and  (d)  the  gas  lubricant  behaves  as  a  perfect  gas  and  compresses  and  ex- 
pands isothermally  within  the  bearing. 

By  introducing  the  boundary  conditions: 

Vr (z=0)  =  0 

Vr(z=h)  =  0 

Ve(z=0)  =  ru> 

Ve(z=h)  =  0 
and  integrating  equations  (1)  and  (2)  the  following  expressions  for  the 
mass  flow  rates  can  be  obtained: 

Qr  =  _  ll.h.  (5) 


Q.=  r[ 


p  u>  h     p  h3    ^ 


12^r*  d9. 


(6) 


Ausman  [4]  introduced  equation  (4)  and  made  the  assumptions  that  (a) 

P<<     may  be  substituted  for   p   ,  ignoring  the  effect  of  compressibil- 

ity  on  the  pressure  gradient  terms  and  (b)  r  '  in  equation  (6)  may  be  re- 

-2 
placed  with  r,   where  r,   is  some  constant  such  that  bR  <   r,  <  R  where 

b  is  the  ratio  of  inner  to  outer  bearing  radii  ■  r^   /R.  The  first  assump- 
tion linearizes  the  resulting  differential  equation  while  the  second  makes 
the  variables  separable.   Further,  introduction  of  two  dimensionless  vari- 
ables, 


r  =  J2 


P- 
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and  substitution  of  equations  (5)  and  (6)  into  equation  (3a)  results 
in  the  linearized,  separable  differential  equation 

A  -  £*# 


where 


Equation  (8)  defines  the  dimensionless  parameter  known  as  bearing  number, 
/\   .    The  form  of  the  solution  is 


P=X(P)T(e)  <9> 


and  results  in  an  infinite  series  solution  which  must  be  solved  indepen- 
dently for  each  flow  region. 

Boundary  conditions  state  that  the  pressure  equals  the  ambient  pres- 
sure at  both  inside  and  outside  bearing  radii,  and  that  the  pressure  and 
flow  rate  are  equal  at  the  step  when  approached  from  either  direction. 
These  latter  conditions  are  met  by  an  iterative  process  which  seeks  a 
value  for  r(   which,  when  substituted  into  the  two  series  solutions,  re- 
sults in  equal  solutions  at  the  step  for  both  regions. 

The  theoretical  load  carrying  capacity  for  a  bearing  is  found  by 
integrating  the  expression  for  p  over  the  bearing  area. 

Ausman  then  optimized  bearing  geometry  by  independently  varying 
step  location,  T/y  ,  and  number  of  pads,  n,  and  solving  for  these  para- 
meters at  maximum  load  capacity.   He  presented  his  results  in  the  form 
of  one  table  and  one  graph  from  which  a  designer  can  determine  the  com- 
plete optimum  geometry  for  a  given  bearing  size  and  expected  operating 
point  (bearing  number,  A).   In  his  optimization  he  assumed  that  the 
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optimum  compression  ratio,  q,  was  2.0  where  <\  =   '/hz   ~  (  hz+  A)/^ 
-   1  +  -^^  .      This,  he  found  from  trial  solutions,  was  a  good  assump- 
tion over  a  wide  range  of  bearing  numbers.  Note  in  Figure  3  (which  is 
described  later)  that  the  solution  for  constant  minimum  clearance  is  maxi- 
mum at  a  clearance  ratio  of  1.0  which  corresponds  to  q  =  2.0. 

Wildmann,  et.    al. ,  in  a  comprehensive  study  [5],  obtained  exten- 
sive experimental  data  at  low  bearing  numbers  which  essentially  con- 
firmed Austrian's  theory  in  the  regions  the  study  undertook  to  investigate. 
Austrian's  results  were  significantly  simplified  by  the  introduction  of  a 
modified  bearing  number 


A'  =  -^fO-b1)  =  A(i-r)  ao) 


This  enabled  him  to  eliminate  the  radius  ratio,  b,  as  an  independent 
parameter  and  to  present  load  capacity  for  optimum  compression  ratio 
(q  =  2)  as  one  curve  plotted  against  bearing  number.  This  curve  is 
shown  as  Figure  2. 

As  a  result  of  experimental  data,  slight  modifications  were  made  to 
Ausman's  theoretical  design  parameters.   New  design  parameters  were  pre- 
sented in  graphical  form  so  that  design  of  a  specific  bearing  could  be 
undertaken  without  the  necessity  of  extrapolating  Ausman's  tabular  re- 
sults. The  first  two  design  charts  are  similar  to  Figures  2  and  3  but 
reflect  Wildmann's  correlation  of  theory.   The  remaining  three  are  plots 
of  optimum  sector  ratio,   T  ■  /)f     ,  and  optimum  sector  arc  length 
(degrees)  and  optimum  number  of  sectors,  n,  all  plotted  vs.  modified 
bearing  number.   The  test  bearing  investigated  in  this  study  is  designed 
using  the  parameters  given  by  Wildmann. 
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One  additional  curve,  computed  from  Ausman's  theory   and  presented 

in  Wildmann's  study,  is  shown  as  Figure  3  because  of  its  usefulness  in 

comparing  experimental  results  with  theory.   This  consists  of  a  plot  of 

minimum  clearance  ratio,  h  3  /  ^      ,  vs. 

,  n  J .      Load 

Load  Ratio  =  -p—~- — i T~i — T~i Z-^ 

Optimized  Load  for  q  =  2 

for  a  constant  step  depth,  A  . 
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3.   Description  of  Apparatus 

Trust  Bearing  Test  Rig 

The  primary  research  tool  used  in  this  experimental  investigation 
is  the  AiResearch  Thrust  Bearing  Test  Rig,  ML  1606,  an  extremely  versatile 
device  designed  to  measure  the  operating  characteristics  of  gas  lubricat- 
ed thrust  bearings  at  very  high  rotation  rates.  Rotating  speed  and  bear- 
ing load  can  be  continuously  varied  and  accurately  measured  while  instru- 
mentation can  be  provided  for  measurement  of  air  film  thickness,  torque, 
and  temperatures  at  selected  points.  Design  variations  permit  testing 
thrust  bearings  up  to  six  inches  in  diameter  and  test  bearings  may  be  of 
any  type  or  configuration. 

Photographs  of  the  Thrust  Bearing  Test  Rig  and  its  associated  instru- 
mentation are  shown  in  Figures  4  and  5.   A  schematic  of  the  arrange- 
ment of  major  components  is  presented  in  Figure  6.   The  rotating  assembly 
consists  of  a  hardened  stainless  steel  shaft  one  inch  in  diameter  with  an 
integral  2%  inch  diameter  thrust  disc  on  the  lower  end  and  a  small  air 
driven  turbine  wheel  on  the  upper  end.  This  shaft  is  supported  vertical- 
ly in  two  AiResearch  ARPF  compliant  journal  bearings.   It  is  supported 
from  beneath  by  the  test  bearing  and  separation  of  the  back  of  the  thrust 
plate  from  the  upper  body  of  the  machine  is  maintained  by  an  externally 
pressuiized  air  thrust  bearing.  This  configuration  is  capable  of  rota- 
tion rates  up  to  105,000  RPM  after  running  through  the  half -frequency 
whirl  region  encountered  in  the  range  between  20,000  and  50,000  RPM. 

The  test  bearing  is  supported  beneath  the  thrust  runner  on  the  top 
face  of  a  heavy  steel  piston  which  moves  vertically  and  through  a  small 
rotation  angle  within  a  cylinder  equipped  with  a  hydrostatic  air  bearing. 

In  the  tests  performed  the  test  bearing  was  recessed  into  the  face 
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of  the  piston  and  was  gimballed  by  placing  it  on  a  3/32  inch  thick  sheet 
of  gasket  rubber.   Restraint  in  torque  was  provided  by  securing  a  .020 
inch  x  \   inch  piece  of  piano  wire  in  a  groove  on  the  piston  face  in  such 
a  way  as  to  allow  the  end  of  the  wire  to  extend  into  the  bottom  of  one 
of  the  air  inlet  grooves  on  the  test  bearing.  A  photograph  of  the  bear- 
ing secured  in  place  is  shown  as  Figure  7. 

Air  Supply  and  Control  System 

A  schematic  of  the  air  supply  system  is  shown  in  Figure  8.   The 
system  utilizes  a  50  HP  Westinghouse  220  VAC,  3  phase,  60  hertz,  Induc- 
tion Motor  to  drive  a  Chicago  Pneumatic,  Model  P2CB  Two  Stage  Air  Compres- 
sor with  intercooler  (max  pressure  =  450  psi).   The  output  air  is  cooled 
by  a  Heliflow  shell  and  tube  heat  exchanger,  Type  9AA14S.  Water  and  oil 
are  removed  by  a  Wright-Austin  Co.,  separator  and  float  trap.  These  are 
followed  by  a  Morse  Diving  Equipment  Co.,  air  filter  and  a  standpipe  con- 
taining a  dessicant  to  further  dry  the  air.  Airflow  to  the  receivers  is 
controlled  by  a  Bailey  Indicator  Controller,  Type  K47,  actuating  a  Bailey 
Diaphram  Valve  Type  \   XSDS-9M.   The  compressed  air  is  stored  in  two 

American  Pipe  &  Steel  Corp. ,  Vertical  Air  Receivers  (each  with  a  volume 

3 
of  117  ft  .  ,  max.  pressure  ■  500  psi).   Air  is  drawn  from  the  top  of 

both  receivers  and  is  distributed  at  outlets  in  a  small  manifold  near  the 
test  setup. 

Figure  9  depicts  schematically  the  air  supply  and  flow  control  ar- 
rangement to  the  test  rig  itself.   Figures  4  and  5  are  photographs  of 
the  actual  apparatus.   Air  is  provided  to  the  driving  turbine  through  a 
one  inch  line  fitted  with  a  regulator  and  globe  valve  for  rough  speed  con- 
trol.  Fine  control  is  achieved  by  a  %  in.  line  equipped  with  a  regulator 
and  needle  valve.   This  line  is  placed  in  parallel  around  the  rough  control 
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line.  A  \   in.  line  with  a  secondary  filter  and  needle  valve  provides 
120  psig  air  to  the  hydrostatic  thrust  bearing.  Air  is  supplied  to  the 
hydrostatic  piston  journal  bearing  at  20  psig  by  a  \   in.  line  with  regula- 
tor and  needle  valve.  A  similar  arrangement  provides  air  to  load  the 
piston  from  beneath,  however  an  atmospheric  bleed  is  placed  in  this  line 
to  prevent  pressure  buildup  due  to  low  flow  rates.   Bourdon  type  pres- 
sure gages  are  located  as  shown  in  the  schematic.   Loading  piston  pres- 
sure is  measured  by  a  14  inch  mercury  U-tube  manometer. 

Electronic  Instrumentation 

Figure  10  schematically  describes  the  electronic  instrumentation 
used  to  measure  rotational  speed,  air  film  thickness  and  temperature. 

Rotational  speed  of  the  test  rig  shaft  is  sensed  by  a  Pace  Engine- 
ering Co.  SU-2  Pickup  which  is  a  very  sensitive  coil  with  a  core  of 
flexible  iron  straps.   The  coil,  when  placed  in  the  vicinity  of  a  rotat- 
ing shaft  of  ferromagnetic  material,  generates  a  small  A.C.  voltage  of 
varying  frequency  as  a  result  of  the  rotating  magnetic  field  due  to  re- 
sidual permanent  magnetism  in  the  shaft.  The  Pace  Pickup  was  placed  as 
shown  in  Figure  4  and  its  output  was  displayed  on  the  y-axis  of  a  Hewlett- 
Packard  Model  130C  Oscilloscope  and  a  Hewlett-Packard  Electronic  Fre- 
quency Meter  Model  500A.   The  output  of  a  Donner  Sine  Wave  Generator 
Model  1202  was  displayed  on  the  X-axis  of  the  oscilloscope  and  its  fre- 
quency was  accurately  measured  by  a  Systron-Donner  Digital  Counter  Model 
1013D.   Since  rotation  rate  is  a  predetermined  figure  for  each  data  run, 
the  rate  i^  preset  on  the  sine  wave  generator  and  the  speed  of  the  rig 
is  adjusted  until  a  stationary  elliptical  lisajous  figure  is  achieved 
on  the  cathode  ray  tube,  indicating  that  the  two  frequencies  are  matched. 
This  results  in  a  speed  measurement  accurate  to  plus  or  minus  one  hertz 
or  +  0.1%  at  60,000  RPM. 
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Air  film  thickness  measurements  were  made  using  Bently  Proximitor 
distance  detectors.   The  calibration  and  accuracy  determination  of  this 
system  is  presented  in  Appendix  A.   During  initial  runs  to  check  out  the 
equipment  it  was  found  that  heating  and  cooling  of  the  distance  probe 
(located  in  the  center  of  the  test  bearing)  was  more  extensive  and  more 
rapid  than  anticipated.   This  change  in  temperature  causes  a  resistance 
change  in  the  sensing  coil  and  hence  changes  the  output  voltage  for  a 
given  distance.   Further,  due  to  low  friction,  rundown  time  (time  neces- 
sary to  stop  the  shaft)  was  in  excess  of  one  minute  and  during  this  time 
the  temperature  and  therefore  the  probe  zero  reading  (thrust  runner  and 
bearing  in  contact),  changed  markedly.  To  correct  this,  two  steps  were 
taken. 

Two  air  nozzles  made  of  \   inch  copper  tubing  were  secured  in  the 
turbine  outlet  throat  to  direct  air  jets  against  the  turbine  blades  op- 
posite to  the  direction  of  rotation.   To  actuate  these,  two  spring  load- 
ed, quick  acting  valves  were  installed.   This  system,  seen  in  Figure  4 
stopped  the  shaft  in  an  average  time  of  15  seconds. 

Two  copper -const ant an  thermocouples  were  installed,  one  adjacent  to 
the  distance  probe  face  and  touching  it,  and  the  other  in  an  ice-water 
bath  as  a  reference  junction.   The  output  was  recorded  on  a  Mosely  Auto- 
graf  Strip  Chart  Recorder  Model  7100B.  The  second  channel  of  the  re- 
corder was  used  as  an  event  marker  to  record  points  at  which  data  was 
taken.  This  arrangement  is  shown  in  Figure  10. 

These  two  corrective  measures  enabled  the  author  to  establish  the 
probe  temperature  difference  between  the  probe  reading  and  the  zero  read- 
ing (thrust  runner  and  bearing  in  contact)  for  each  data  point.   The 
significance  of  this  data  will  be  shown  in  Appendix  C  -  Experimental 
Uncertainty  Intervals. 
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Test  Bearings 

It  was  originally  intended  to  test  several  thrust  bearings  of  dif- 
ferent optimum  configurations  but  time  and  funding  limitations  restricted 
the  tests  to  one  bearing.  The  bearing  was  chosen  to  correspond  to  the 
point  in  Austrian's  Table  I  [4]  at  b  -  .20,  .A  -  A  ^  =  160.  The  design 
parameters  however,  were  taken  from  Wildmann  [5].  These  parameters  are: 

b  =  .20 
Al-A*(l-kl)«  153-7 

n  «  3 

V  =  118° 

*    =  16° 
The  minimum  film  thickness  was  chosen  to  be  250.  microinches  so,  for 
q  ■  2.0,  the  optimum  step  depth,   A   ,  is  250.  microinches.   The  out- 
side radius,  R,  was  chosen  as  1.125  inches.  The  final  design  based  on 
these  parameters  is  shown  in  Figure  11.   Note  the  air  inlet  groove  which 
is  two  degrees  wide  at  the  average  radius.   The  material  chosen  for  the 
bearing  was  "Graphmo"  tool  steel,  a  graphite-molybdenum  steel  noted  for 
its  wear  resistance  and  low  sliding  friction.  The  bearing  was  manufactur- 
ed to  specifications  by  B  &  W  Precision,  Inc.,  Glendale,  California. 

Inspection  of  the  bearing  by  optical  flat  before  running  revealed 
a  flatness  of  +  5.0  microinches  across  the  bearing  face.   Profilometer 
readings  by  the  manufacturer  showed  a  surface  roughness  of  5  microinches 
(rms) .   The  step  depth  was  measured  on  a  Johansson  Mikrokator  Comparator 
and  varied  between  230  and  270  microinches  with  an  average  of  250.   These 
readings  were  taken  numerous  times  across  all  three  steps.   The  surface 
of  the  recessed  area  was  somewhat  rougher  than  the  raised  portion  but  was 
not  measured. 
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4.    Experimental  Procedure 

Data  was  obtained  in  such  a  manner  that  load  vs.  air  film  thickness 
could  be  plotted  for  a  single  constant  speed.   This  was  repeated  for  a 
number  of  different  speeds  to  obtain  a  family  of  such  curves.   No  data 
was  taken  in  the  region  where  half-frequency  whirl  was  encountered  to 
insure  that  the  data  was  not  affected  by  this  instability. 

Test  Rig  Operation 

All  runs  were  made  with  the  hydrostatic  thrust  air  set  at  115  to 
120  psig  and  the  hydrostatic  piston  journal  air  adjusted  to  20  psig. 
All  starts  and  stops  were  performed  with  a  test  bearing  load  of  3.65 
pounds  (3.0  inches  Hg) .   Prior  to  starting,  the  speed  for  the  run  (re- 
volutions per  second)  was  preset  on  the  digital  counter  by  adjusting  the 
sine  wave  generator  to  the  same  frequency  (hertz) .   The  turbine  was 
then  started  by  quickly  opening  the  rough  control  valve  enough  to  bring 
the  downstream  pressure  to  about  5  psig.   The  bearing  noise  was  monitor- 
ed continuously  for  whirl  or  rubbing  by  listening  with  an  automotive  type 
stethoscope.   Rotational  speed  was  adjusted  while  watching  the  lisajous 
figure  on  the  oscilloscope  and  reading  the  frequency  meter  until  the 
speed  was  slightly  greater  than  the  preset  rate  and  steady.  The  tempera- 
ture recorder  chart  was  then  started  and  the  bearing  load  slowly  increas- 
ed, taking  care  that  the  temperature  did  not  increase  at  a  rapid  rate  and 
cause  excessive  thermal  distortion.  As  the  bearing  was  loaded,  the  tur- 
bine air  flow  was  increased  to  maintain  speed. 

When  the  load  for  the  particular  data  point  was  reached,  the  rig 
was  allowed  to  run  until  the  temperature  reached  a  steady  state  or  was 
increasing  very  slowly.   Then  the  fine  speed  control  flow  valve  was  used 
to  control  turbine  speed  until  the  oscilloscope  trace  showed  an  exact 
frequency  match.   At  this  point,  the  distance  probe  output  was  read  and 
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noted,  the  turbine  air  was  secured,  the  braking  jets  were  actuated,  and 
the  bearing  load  was  reduced  to  3.65  pounds.  As  the  probe  output  was 
read,  the  event  was  marked  on  the  strip  chart  recorder.  At  the  instant 
that  the  shaft  stopped  and  the  thrust  disc  and  bearing  came  in  contact, 
the  distance  probe  output  was  again  read  and  the  event  marked. 

This  procedure  was  repeated  for  each  data  point  obtained,  and  each 
data  point  was  repeated  several  times  to  insure  correctness  and  to  ob- 
tain a  statistical  average.  The  data  obtained  are  presented  in  Tables 
1  through  4. 

Bearing  Failure 

On  the  last  run  for  data,  a  failure  occurred  and  the  thrust  plate 
contacted  the  bearing  face  with  subsequent  rapid  loss  of  air  film  and 
wiping  of  the  bearing  surfaces.  The  failure  was  initiated  at  a  speed  of 
72,000  RPM  and  a  bearing  load  of  14.7  pounds.  Just  prior  to  failure  a 
vibration  noise  began  and  increased  in  amplitude  even  though  the  load 
was  reduced  sharply.  The  center  of  both  the  bearing  and  thrust  disc 
were  wiped  to  a  radius  of  approximately  .35  inches  which  would  seem  to 
indicate  failure  due  to  thermal  distortion.   This  is  felt  not  to  be  the 
case,  however,  since  neither  temperature  nor  temperature  rate  were  ex- 
cessive at  the  time,  and  the  load,  while  high,  was  no  higher  than  for 
several  other  runs  when  the  air  film  thickness  was  even  less.   Rather 
it  is  felt  that  the  speed  used  was  a  critical  vibration  frequency  for 
some  component,  most  likely  the  rubber  sheet  used  to  gimbal  the  test 
bearing.  The  appearance  of  the  surfaces  after  failure  merely  indicates 
that  slight  thermal  distortion  was  present  at  the  time  of  failure. 

Data  Reduction  Calculations 

The  data  reduction  calculations  used  to  obtain  the  results  presented 
in  the  next  section  are  described  in  Appendix  B. 
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5.    Results 

The  numerical  results  calculated  from  the  data  obtained  are  present- 
ed in  Table  5. 

Figures  12  through  15  are  plots  of  bearing  load  vs.  air  film  thick- 
ness for  each  speed  tested.  All  data  points  are  plotted  and  extreme 
uncertainty  limits  are  shown.   Figure  16  is  a  composite  plot  of  the  same 
data  using  only  the  statistical  averages  for  air  film  thickness. 

Figure  17  is  a  graph  of  dimensionless  load  parameter,  L,  vs.  the 
modified  bearing  number,   -A-jfi  f°r  a  compression  ratio,  q  ■  2.0. 

Perhaps  the  most  useful  form  of  presenting  the  results  is  seen  in 
Figure  18.   The  ratio  of  actual  load  to  optimized  load  at  q  =  2.0  is 
shown  on  the  ordinate  logarithmically  while  the  clearance  ratio, ^  ■  y^ 
is  plotted  on  the  abcissa.   By  the  optimized  load  it  is  meant  the  opti- 
mized load  from  the  theoretical  curve  in  Figure  2  determined  for  the 
same  bearing  number  as  the  particular  data  point  under  consideration.  A 
clearance  ratio  of  1.0  corresponds  to  q  ■  2.0.  All  of  the  data  points 
from  Figure  17  are  shown  in  Figure  19. 
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6.   Comparison  of  Results  with  Theory 

In  Figure  19,  the  results  shown  in  Figure  17  are  superimposed  on  a 
reproduction  of  Wildmann's  Design  Chart  A  [5],  which  predicts  load  carry- 
ing capacity  at  q  ■  2.0  for  a  given  bearing  number,  TV   .  Note  that 
the  results  take  roughly  the  same  form  as  the  predicted  curve  but  the 
load  capacity  is  signifigantly  lower. 

In  Figure  20,  the  curve  from  Figure  18  is  replotted  on  a  graph  which 
shows  both  Ausman's  theoretical  prediction  and  Wildmann's  experimental 
results.  Note  that  the  results  of  this  study  have  not  been  normalized 
to  a  load  ratio  of  1.0  as  was  done  in  the  study  compiled  by  Wildmann. 

Two  factors  are  immediately  apparent.  The  first  is  that  the  slope 
of  the  curve  obtained  in  the  current  study  is  essentially  the  same  as 
Wildmann's  correlation  curve.   Secondly,  the  new  curve  predicts  a  load 
capacity  that  is  a  constant  fraction,  2/5  or  0.4  of  that  predicted  by 
the  design  curves  [f5]  based  on  Ausman's  work. 
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7.   Discussion  of  Results 

It  is  felt  that  the  careful  calibration  of  the  probes  used  to  meas- 
ure air  film  thickness,  plus  the  temperature  monitoring  employed,  re- 
sulted in  reliable  air  film  thickness  measurements.   Further,  the  repeti- 
tion of  data  points  and  use  of  average  values  reduced  random  error.   This 
opinion  is  borne  out  by  the  smoothness  of  the  load  vs.  film  thickness 
curves  and  the  consistency  of  the  load  ratio  vs.  clearance  ratio  curve. 

There  are  two  sources  of  possible  major  error.   The  first  is  due  to 
the  near  impossibility  of  constructing  the  perfect  stepped  sector  thrust 
bearing  upon  which  the  theory  is  based.   It  was  noted  that  the  variation 
in  the  step  depth  was  plus  and  minus  20  microinches  or  127.  of  the  speci- 
fied depth.   This  plus  the  rougher  finish  on  the  recessed  area  undoubted- 
ly affected  the  flow  characteristics  of  the  lubricant.  The  second  seri- 
ous condition  is  thermal  distortion.   Even  though  attempts  were  made  to 
allow  temperature  to  reach  a  steady  state,  distortion  was  surely  present 
as  shown  by  the  nature  of  the  bearing  failure  which  occurred.  Thermal 
distortion  would  result  in  a  nonuniform  air  film  thickness;  the  film 
would  increase  in  thickness  as  the  radius  at  which  it  is  measured  is  in- 
creased.  It  is  not  felt,  however,  that  either  of  these  conditions  are 
responsible  for  the  wide  disparity  between  theory  and  experiment  with 
regard  to  load  capacity,  particularly  since  the  data  is  otherwise  con- 
sistent. 
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8.   Conclusions 

From  the  comparison  of  results  with  theory,  it  is  concluded  that: 
(a)  Wildmann's  correlation  of  the  variation  of  load  ratio  with  clearance 
ratio  is  correct  at  high  bearing  numbers  for  the  bearing  tested;  (b) 
Ausman's  theory  correctly  predicts  the  shape  of  the  load  vs.  bearing 
number  curve  in  the  high  bearing  number  region,  and  (c)  the  prediction 
of  the  magnitude  of  the  load  capacity  is  in  error  in  this  region. 

Ausman,  in  the  preparation  of  his  design  table,  performed  separate 
optimizations  for  bearing  geometry  but  extrapolated  load  capacity  from 
curves  computed  for  a  central  point  in  the  table.   He  warned  his  readers 
that  the  accuracy  of  this  extrapolation  probably  decreased  with  distance 
from  the  center  of  the  table.  The  results  of  this  test,  of  a  bearing 
purposely  chosen  from  an  extreme  corner  of  his  table,  confirm  that  Ausman's 
intuition  was  correct. 
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9.    Recommendations  for  Further  Study 

Now  that  it  is  established  that  there  is  a  limit  to  the  range  of 
accuracy  of  Austrian's  design  parameters,  additional  bearings,  optimized 
by  the  procedure  used  in  this  study,  should  be  tested  until  the  range 
of  accuracy  is  established. 

A  computer  aided  analysis  of  the  load  capacity  equations  should  be 
performed  over  the  full  range  of  bearing  numbers  using  the  modified 
bearing  number  introduced  by  Wildmann.   The  results  should  give  some 
insight  into  the  true  nature  of  the  relationship  and  could  result  in  a 
revision  of  the  optimized  design  parameters. 

Other  types  of  thrust  bearings,  such  as  the  spiral  groove  or  her- 
ringbone patterns,  await  testing  at  high  bearing  numbers.   The  AiResearch 
Thrust  Bearing  Test  Rig  is  ideal  for  this  purpose. 
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TABLE  1 
Summary  of  Data  --  54,000  RPM 


RUN  NO. 

14-1 
14-2 
14-3 
14-4 
14-5 
14-6 
15-1 
15-2 
15-3 
15-4 
15-5 
15-6 
16-1 
16-2 
16-3 
16-4 
16-5 
16-6 
17-1 
17-2 
17-3 
17-4 
17-5 
17-6 


LOAD 

h* 

(pounds) 

(^in.) 

6.59 

363 

6.59 

366 

6.59 

371 

6.59 

363 

6.59 

366 

6.59 

370 

8.62 

318 

8.62 

330 

8.62 

337 

8.62 

320 

8.62 

325 

8.62 

320 

10.65 

287 

10.65 

287 

10.65 

287 

10.65 

284 

10.65 

289 

10.65 

294 

12.67 

254 

12.67 

260 

12.67 

256 

12.67 

258 

12.67 

254 

12.67 

253 
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TABLE  1  (Cont) 

RUN  NO.                  LOAD  h^ 

(pounds)  ( ^  in.) 

18-1                    13.69  246 

18-2                    13.69  246 

18-3                    13.69  244 

18-4                    13.69  251 

18-5                    13.69  243 

18-6                    13.69  250 
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TABLE  2 
Summary  of  Data  --  60,000  RPM 
RUN  NO.  LOAD 


(pounds) 


2 
(^in.) 


1-1  no  data 

2-1  4.06  442 

2-2  4.06  452 

2-3  4.06  455 

2-4  4.06  466 

2-5  4.06  456 

2-6  4.06  455 

3-1  6.59  360 

3-2  6.59  363 

3-3  6.59  370 

3-4  6.59  380 

3-5  6.59  385 

3-6  6.59  394 

4-1  8.62  320 

4-2  8.62  313 

4-3  8.62  330 

4-4  8.62  333 

4-5  8.62  337 

4-6  8.62  328 

5-1  10.65  298 

5-2  10.65  294 

5-3  10.65  294 

5-4  10.65  296 

5-5  10.65  299 
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RUN  NO. 

5-6 
6-1 
6-2 
6-3 
6-4 
6-5 
6-6 
7-1 
7-2 
7-3 
7-4 
7-5 
7-6 
7-7 
8-1 
8-2 
8-3 
8-4 
8-5 
8-6 


TABLE  2  (Cont) 

LOAD 
(pounds) 

h2 
(<\in.) 

10.65 

289 

12.67 

256 

12.67 

257 

12.67 

256 

12.67 

268 

12.67 

261 

12.67 

260 

13.69 

243 

13.69 

251 

13.69 

253 

no  data 

13.69 

253 

13.69 

251 

13.69 

253 

14.19 

240 

14.19 

250 

14.19 

241 

14.19 

240 

14.19 

243 

14.19 

240 
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RUN  NO. 

9-1 

9-2 

9-3 

9-4 

9-5 

9-6 
10-1 
10-2 
10-3 
10-4 
10-5 
10-6 
11-1 
11-2 
11-3 
11-4 
11-5 
11-6 
12-1 
12-2 
12-3 
12-4 
12-5 


TABLE  3 

Summary  of  Data  --  66,420  RPM 

LOAD 
(pounds) 

h2 
(>\   in.) 

6.59 

376 

6.59 

380 

6.59 

385 

6.59 

385 

6.59 

377 

6.59 

374 

8.62 

332 

8.62 

347 

8.62 

337 

8.62 

339 

8.62 

336 

8.62 

337 

10.65 

299 

10.65 

305 

10.65 

312 

10.65 

295 

10.65 

295 

10.65 

298 

12.67 

263 

12.67 

274 

12.67 

264 

12.67 

274 

12.67 

268 
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TABLE  3  (Cont) 
RUN  NO.  LOAD  h 


(pounds) 


2 
0\in.) 


13-1  14.70  250 

13-2  14.70  247 

13-3  14.70  240 

13-4  14.70  247 

13-5  14.70  243 

13-6  14.70  246 
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TABLE  4 

Summary  of  Data  —  72,000  RPM 

LOAD 
(pounds) 

h2 
(  H  in.) 

8.62 

337 

8.62 

353 

8.62 

346 

10.65 

301 

10.65 

301 

10.65 

301 

12.67 

265 

12.67 

267 

12.67 

277 

13.69 

258 

13.69 

260 

14.70 

253 

RUN  NO. 

19-1 
19-2 
19-3 
20-1 
20-2 
20-3 
21-1 
21-2 
21-3 
22-1 
22-2 
23-1 
23-2  14.70  Bearing  failure. 
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SECTION  A -A 
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NOT  SHOWN  TO  SCALE  —  CLEARAHCtS  GREATLY  SXMDGER4TVD 


FIGURE  1.    TYPICAL  STEP  BEARING  GEOMETRY 
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FIGURE  6.    MAJOR  COMPONENTS  OF  THRUST  BEARING  TEST  RIG 
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FIGURE  12.    BEARING  LOAD  vs.  AIR  FILM  THICKNESS  -  54,000  RPM 
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FIGURE    13.        BEARING   LOAD   vs.   AIR  FILM  THICKNESS   -   60,000  RPM 
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FIGURE  14.    BEARING  LOAD  vs.  AIR  FILM  THICKNESS  -  66,420  RPM 
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FIGURE  15.    BEARING  LOAD  vs.  AIR  FILM  THICKNESS  -  72,000  RPM 
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FIGURE  16.   BEARING  LOAD  vs.  AIR  FILM  THICKNESS  AT  VARIOUS  SPEEDS 
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APPENDIX  A 
DISTANCE  PROBE  CALIBRATION 
Selection  of  Measuring  Systems 

The  thickness  of  the  air  film  between  the  test  bearing  and  its 
thrust  runner  was  expected  to  vary  between  150  to  400  microinches.   It 
was  felt  that  significant  data  could  fee  obtained  with  a  measuring  techni- 
que which  could  discriminate  within  olus  or  minus  10  microinches  during 
actual  operation.   Based  on  descriptions  of  apparatus  found  in  the  liter- 
ature, two  types  of  sensors  should  be  considered  ---  a  capacitance  probe 
which  would  use  the  air  film  to  be  measured  as  the  dielectric  between 
two  parallel  plates,  and  the  magnetic ,  or  reluctance,  type  which  senses 
the  proximity  of  a  conductor  by  the  extent  to  which  it  disturbs  a  mag- 
netic field  surrounding  a  sensing  coil. 

Although  it  is  capable  of  high  accuracy,  the  capacitance  type  was 
eliminated  because  it  depends  on  the  dielectric  constant  of  the  air  film. 
It  was  felt  that  changes  in  temperature  and  moisture  content  would  un- 
predictably alter  the  dielectric  properties  of  the  air  film. 

Bently  Proximitors  manufactured  by  the  Bent ly -Nevada  Corp. were  chosen 
because  of  their  high  potential  accuracy,  convenient  mounting  arrange- 
ments, and  low  cost.   Since  they  are  of  the  magnetic  reluctance  tvpe, 
their  measurements  depend  upon  the  nature  c  *  the  conducting  surface  being 
sensed.  Variations  in  the  properties  of  nonconducting  materials  be- 
tween the  probe  and  the  measured  surface  have  no  measurable  effect  on 
the  accuracy  obtained.   Temperature  of  the  probe  itself  was  found  to 
have  an  effect,  however,  even  though  silver-wound  coils,  advertised  as 
the  least  temperature  sensitive,  were  used. 
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Enhancement  of  Accuracy 

The  selected  probes  were  intended  to  measure  a  range  of  0  to  .010 
inches  and  to  resolve  to  +  .0001  inches.  To  attain  the  necessary  ac- 
curacy (+  10  microinches) ,  the  resolution  had  to  be  improved  tenfold. 

According  to  the  manufacturer,  the  accuracy  and  resolution  is  direct- 
ly proportional  to  the  stability  of  the  input  voltage.  Accordingly,  an 
ultra-stable  power  supply  meeting  the  voltage  and  current  requirements 
was  purchased  from  Systems  Research  Corporation.  Measured  ripple  and 
noise  of  their  Model  3566  Power  Supply  was  found  to  be  50  to  80  micro- 
volts peak  to  peak  and  the  output  varied  no  more  than  1.0  millivolts  in 
a  24  hour  test. 

To  read  the  output  voltage  accurately,  a  digital  voltmeter,  Electro 
Instruments  Co.  Model  DVX-500,  capable  of  measuring  to  the  nearest  0.1 
millivolts  was  obtained. 

Probe  Calibration  Methods 

Photographs  of  the  probe  calibration  rig  used  are  shown  in  Figures 
Al  and  A2.   The  apparatus  consists  of  a  Johannson  Mikrokator  Model  S 
510E-7  Comparator  altered  for  use  as  a  calibrator.   The  Mikrokator  dial 
indicator  seen  in  Figure  A2  is  accurate  to  the  nearest  one  millionth  of 
an  inch  and  is  thus  capable  of  establishing  the  actual  resolution  of  the 
distance  probes.   The  comparator  was  altered  by  providing  a  rigid  clamp 
fixture  to  hold  a  distance  probe  over  the  movable  table,  and  by  making 
arrangements  to  adjust  the  position  of  the  dial  indicator  head  by  means 
of  a  fine  adjustment  mechanism  salvaged  from  a  microscope. 

The  surface  to  which  the  probes  were  calibrated  was  the  surface  of 
a  disc  made  indentical  to  the  thrust  runner  actually  used  in  the  thrust 
bearing  test  rig.   This  was  done  by  cutting  the  disc  from  the  same  billet 
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of  steel  and  having  it  follow  the  shaft  with  thrust  runner  through  the 
same  machining,  heat  treating,  grinding  and  lapping  processes. 

Provisions  were  made  for  checking  the  effects  of  temperature  on  the 
distance  measurements.   The  ambient  temperature  was  accurately  recorded 
before,  during,  and  after  each  calibration.   The  calibration  rig  was  al- 
lowed to  "soak"  for  at  least  three  hours  at  any  given  temperature  to  elim- 
inate any  transient  expansion  or  contraction. 

A  temperature  sensitivity  test  was  performed  to  insure  that  tempera- 
ture changes  in  readout  were  not  due  to  mere  distance  changes  as  a  result 
of  expansion  or  contraction  of  clamping  fixtures  of  the  calibration  rig 
itself.   A  probe  was  clamped  in  a  vertical  position  with  its  sensing  sur- 
face up.  The  test  surface  disc  was  allowed  to  rest  above  the  probe  iace 
with  a  separation  of  .015  inches  maintained  by  a  thin  sheet  of  Pyrex 
glass.  The  coefficient  of  thermal  expansion,  °*  ,  of  Pyrex  glass  is 
2.0x10   -: — 10    So  for  a  temperature  change  of  100F,  for  instance,  the 
glass  thickness  ,a  ,  will  change  by 

A  ^  _  w  •  [  i\  i  j 

-  (2.0xl0"6)(.015)(100) 

=  3.0  microinches 
This  arrangement  was  then  subjected  to  heated  air  blown  from  a  portable 
heat  gun  while  the  output  reading  was  watched  for  deviation. 

Results 

In  the  temperature  sensitivity  test,  raising  the  temperature  ap- 
proximately 50F  caused  a  change  in  reading  of  about  100  millivolts.   This 
represents  a  relative  distance  change  of  about  140  microinches  and  con- 
clusively confirmed  suspicions  that  the  probes  themselves  are  affected  by 
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temperature. 

Calibration  curve  data  are  presented  in  Tables  A1-A3.   The  resulting 
calibration  curves  are  presented  in  Figure  A3.   Note  that  data  was  taken 
only  in  a  small  linear  region  of  the  probe's  effective  range  since  larger 
displacement  measurements  were  not  anticipated.  The  slope  of  the  curves 

was  the  same  in  every  case: 

.,      i  /iio  microinches 

Slope  ■  1.4118     . r- — 

millivolt 

Discussion 

From  Figure  A3  it  is  immediately  apparent  that  although  temperature 
changes  affect  the  position  of  the  curve,  only  the  zero,  or  reference, 
point  is  affected.  The  slope  of  the  curve  remains  constant  over  the  ex- 
pected range  of  measurements.  Accordingly,  the  experimental  procedure 
during  actual  runs  included  the  determination  of  a  zero  point  (bearing 
and  thrust  runner  in  contact)  immediately  after  each  run  while  the  rig 
and  probe  were  still  at  operating  temperature.  Measurements  were  then 
taken  from  this  zero  point  using  the  constant  slope  of  the  calibration 
curves. 

To  determine  data  scatter  and  accuracy,  the  data  was  originally 
plotted  on  20x24  inch  graph  paper  with  10x10  squares  to  the  half -inch. 
The  mean  deviation  from  the  straight  line  plots  was  0.96  microinches. 
Thus  by  accuracy  and  stability  enhancements,  together  with  careful  cali- 
bration, the  resulting  data  may  be  considered  reliable  within  +1.92  micro- 
inches,  or  twice  the  mean  deviation.   This  is  better  than  five  times  the 
necessary  accuracy  and  fifty  times  better  than  the  manufacturer's  intended 
accuracy. 
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TABLE  Al 
DISTANCE  PROBE  CALIBRATION  DATA 
68,6°F 
Run  No,  Relative  Distance  Probe  Output 

(in.)  (rav) 

1  400  6069 

379  6035 

337  6000 

296  5968 

268  5945 

204  5900 

170  5875 

121  5838 

66  5800 

3  5756 

2  392  6056 

359  6031 

320  6005 

271  5968 

216  5933 

168  5896 

111  5857 

63  5820 

9  577? 
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TA3LE  A2 
DISTANCE  PROBE  CALIBRATION  DATA 
75.4°F 

Run  No.                  Relative  Distance  Probe  Output 

(in.)  (mv) 

3  398  5943 

350  5907 

302  5874 

251  5836 

193  5795 

139  5759 

95  5724 

52  5695 

4  398  5941 

344  5898 

310  5872 

229  5815 

176  5777 

114  5731 

62  5689 

0  5648 
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TABLE  A3 
DISTANCE  PROBE  CALIBRATION  DATA 
49.4°F 
Run  No,  Relative  Distance  Probe  Output 

(in.)  (mv) 

5  396  6242 

366  6219 

307  6175 

254  6134 

234  6123 

177  6082 

217  6037 

68  6005 

0  5956 

6  395  6207 

360  6186 

310  6154 

251  6117 

207  6088 

158  6054 

123  6028 

87  6004 

51  5981 

18  5956 

0  5942 
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APPENDIX  B 
DATA  REDUCTION  CALCULATIONS 
The  results  of  measurements  shown  in  the  Summary  of  Data  (Tables 
1  through  4;  were  calculated  from  the  raw  data.   Examples  of  these  calcu- 
lations follow. 

Since  the  speed  measurement  is  indicated  by  a  frequency  (hertz)  and 
1  hertz  ■  1  revolution  per  second 
speed  in  revolutions  per  minute  is  computed  by  multiplying  this  value  by 
60.   For  example: 

N  -  900  hz.  x  60  =  54,000  rpm 

The  bearing  load  is  indicated  by  the  height  of  a  column  of  mercury 
in  inches.  The  manometer  reading  is  converted  to  psi  and  multiplied  by 
the  piston  area.  The  weight  of  the  piston  and  test  bearing  assembly  is 
subtracted  from  the  result  to  obtain  the  actual  bearing  load: 

2.036  in.  Hg  =  1.0  psi 
Piston  Area  =  10.317  in. 
Weight  of  Piston  Assembly  =  11.65  lbs. 
Therefore,  for  a  reading  of  4.0  in.  Hg, 

n  -  4.0  x  10.317   ..  ,_   Q  ,_  ., 

W  ■  .  _-, -  11.65  =  8.62  lbs. 

Z.OJo 

Air  film  thickness  is  determined  from  the  difference  (in  millivolts) 

between  the  voltmeter  reading  at  a  particular  speed  and  load  and  the 

reading  immediately  thereafter  taken  with  the  bearing  and  thrust  disc  in 

contact.   This  difference  is  multiplied  by  the  constant  slope  of  the 

probe  calibration  curve  (Appendix  A)  to  obtain  air  film  thickness.   From 

Run  No.  14-1, 
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h  =  (5864  -  5607)(1.4118) 

=  257  x  1.4118  =  363  microinches 
Subsequent  analysis  of  the  results  is  performed  using  the  dimension- 
less  parameters,        ,  L,   -   ,  and  Load  Ratio  (L/L     ).    The  calcu' 
lation  of  these  parameters  is  described  in  the  following  paragraphs. 
The  bearing  number,        ,  is 


where  constant  values  are 


'  =  28  x  10"10  reyns 


R  =  1.125  in, 

=  25.0  x  10"!'  in. 


b  =  0.2 


and 


For  Run  No.  6, 


and 


p  =  barometric  pressure  (in.  Hg)  x  ,4912  ~i~"7\ 

Cn,  n 


=  1000  (2  )  =  6283  rad./sec. 
p  =  30.31  x  .4912  =  14,89  psia 

A'  _  (6)  (28  x  10"10)(6283)(1.1252>   «    22 
(14.89) (25  x  UT*)1  v    '   ; 


=  137.8 
The  load  parameter,  L,  is  given  by 

p,  A 
where 


A  =   (R2    -   rt2)7rx  (354/360) 

=   (1.1252   -    .2252)    r   (354/360) 


=  3.755   in. 
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The  fraction  354/360  accounts  for  the  fact  that  the  air  inlet  grooves 
are  not  part  of  the  load  carrying  area.   Again  for  Run  No.  6, 

t  _     12.67 


(14.89)(3.755) 
=  0.227 
The  clearance  ratio,   c<  =  h./A    »  is  computed  from  the  average 
film  thickness  for  a  given  run,   Using  Run  No.  6, 


and 


h  (avg)  =  259.7  microinches 


•<  -  SH  ■  *•« 


The  optimized  load  for  q  =  2  is  found  by  consulting  Figure  3  for 
the  particular  bearing  number,  -A.^  .   For  example  if 

Al=  137.8 
L  =  .628 

and 

.  227 
Load  Ratio  ■  Zoo"   =  »361 
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APPENDIX  C 
EXPERIMENTAL  UNCERTAINTY  INTERVALS 
Reference | 9  gives  the  uncertainty  interval  for  a  result  as 


■ 


where  w   ,  w   ,...,  w   are  the  uncertainty  intervals  of  the  variables 
making  up  the  result,   R  is  the  expression  for  the  result,  and  v    ,   V 
...,  v    are  the  variables  making  up  the  result. 
The  expression  for  bearing  load  is 


■ 


1 


'! 


where 


P 

A 

i 

W 


■  load  pressure  (in.  Hg) 

=  piston  area  (in.) 

=  weight  of  piston  &  test  bearing  assembly  (lbm.) 


The  uncertainty  interval  for  bearing  load  is 


. 


Dividing  through  by  W  to  nondimensionalize  the  equation  results  in 


j 


Lx 


w. 


t- 


-, 


K 


+  i 


I   ^w* 


.  1  /a 


v  H  U 


For 


P:    =  4-4  1  *05  in-  H§ 
Af   =  10.317  +  .001  in. 

=  11.65  +  .02  lbm 


(20  to  1) 
(20  to  1) 
(20  to  1) 


ii- 


-  .01544 
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and 

w    =  (10.65) (.01544)  =  +  0.164  lbf. 


v7 


Air  film  thickness  is  computed  in  the  following  way: 

h2  =  (V1  "  V2)S 
where 

V-  =  probe  output  voltage  at  speed  and  load  (mv) 

V  =  probe  output  voltage  at  zero  point  (mv) 

S  ■  slope  of  probe  calibration  curve  (m   in./mv) 

One  additional  factor  must  be  taken  into  account.   Inspection  of  the  probe 

calibration  curves  shows  an  average  temperature  sensitivity  of  14  mv/°F 

at  a  constant  distance  between  probe  and  measured  surface.   Examination 

of  all  temperature  data  taken  shows  that  in  only  one  out  of  ten  readings 

did  the  probe  temperature  differ  at  speed  and  load  from  the  zero  point 

by  more  than  1.25  °F.   An  actual  temperature  trace  from  the  strip  chart 

recorder  is  shown  as  Figure  C-l.   Defining  an  additional  error, 

w    =  uncertainty  interval  due  to  a  difference  in  prcfee 
temperatures 

=  (14.0)(1.25) 

=  +  17.5  mv 
It  is  felt  by  the  writer  that  repetition  of  each  data  point  an  average  of 
six  times  halved  this  uncertainty  interval  at  the  very  least.   The  inter- 
val used  is 

w    =  +  8.75  mv      (10  to  1) 
Therefore  for  Run  No.  5-2  using  the  previous  method, 


\j-1  'Z2- 


where 
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V,    =  5788  +  2  mv    (20  to  I) 


and 


V^   =  5580  +  2  mv    (20  to  1) 

S    =  1.4118  +  .005,,in./mv  (20  to  1) 

w,.   ■  294  (.0422) 

='12.41  microinches 
The  expression  for  bearing  number  is 

V-  I- 

;  ,  A'" 


(!-!') 


and  by  similar  means 

where  for  Run  No.  5-2, 

A[     -  28xl0"10  +  lxlO"10  reyns  (20  to  1) 

A     =  6283  +  6.283  rad./sec.  (20  to  1) 

p,  ■  30.31  +  .01  in.  Hg  (20  to  1) 

A  =  250  +  10  microinches  (20  to  1) 

wA   =  137. 8(. 0669) 
=  +  9.22 
Similarly  for  Run  No.  5  (averaged  points) , 

/a 


and 


ft. 
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VL 


A 


ji 


7~/T 


2<7<r 


ft 


w^   -  1.18  (.0581) 
=  +  .069 
The  foregoing  are  examples  of  the  determination  of  uncertainty 
intervals  for  each  data  point.   These  intervals  are  represented  on  all 
graphical  presentations  of  data  as  short  lines  through  the  data  points  in 
the  following  manner: 
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FIGURE  CI.    PROBE  TEMPERATURE  AS  A  FUNCTION  OF  TIME 
FOR  RUN  NO.  17-4 
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